Effect of Light Fermions on the Confinement Transition in QCD-like Theories 
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Dependence of the confinement transition parameters on the fermion content provides information 
on the mechanism of confinement. Recent progress in lattice gauge theories has allowed to study it 
for light flavor number Nf ~ O(10) and found this transition to shift toward significantly stronger 
coupling. We propose an explanation for that: light fermions can occupy the chromo-magnetic 
monopoles, via zero modes, making them "distinguishable" and unsuitable for Bose-Einstein Con- 
densation. Such dilution of unoccupied monopoles is compensated by stronger coupling that makes 
them lighter and more numerous. We also suggest that flavor-carrying quark-monopole objects 
account for the density beyond quark Fermi sphere seen in cold dense phase of N c — 2 lattice QCD. 



1. The color confinement of Quantum Chromodynam- 
ics (QCD) remains to be one of the most outstanding 
puzzles of the Standard Model, in spite of intense studies 
via lattice QCD simulations |T| . By virtue of the asymp- 
totic freedom, at high temperature the effective coupling 
is weak and the matter is in a deconfined phase known as 
the quark-gluon plasma, or QGP 2\. When T is lowered, 
as has happened after the Big Bang or in heavy ion col- 
lisions, the effective coupling grows. At certain critical 
temperature T c , a transition into the confined hadronic 
world occurs [31 H]. Physics had many examples in which 
the dependence of T c on some parameters had offered 
crucial insights: e.g. when the T c of superconductivity 
had shown an isotope mass dependence. Our strategy is 
similar: examining how this transition depends on the 
fermion representation and flavor number Nf may lead 
to insights about the mechanism of confinement. 

Our phcnomcnological input comes from the lattice 
studies. A (very incomplete) list of those ranges from 
the well studied region of Nf = 2,3 (see e.g. E] to the 
recent extension toward Nf — 8 (7] [8] and even Nf = 12 
[S], which has attracted special attention in connection 
with the search for a conformal regime, see e.g. [10] . 
There are also studies with adjoint [TTHT5] and tensor- 
symmetric quarks [14] in similar regime, see e.g. reviews 
[T51 HB] • Starting from the pure gauge theory (Nf = 0) 
and increasing Nf, one finds the monotonous and persis- 
tent shift toward the stronger coupling at T c . (The value 
of T c itself is usually expressed via units customary in 
lattice community, which fixes the T = string tension 
to the same real-world value. We will not use such units 
as they confuse the comparison across different QCD-like 
theories.) The absolute magnitude of the gauge coupling 
constant f3 c at the transition temperature T c , instead, 
bears more direct information. For doing so, we evolve 
the critical lattice coupling /3l — 2Nc/g 2 at the lattice 
scale a (by 2-loop running) to the scale 1/T C = N T a : 

[■—) e 2W^Pc- N 2 f-^-J e 2N c „ fe (1) 
where 6q, 6i are the usual /3-function coefficients. In FigJT] 



we've collected such j3 c values for theories with varied Nt 
from various lattice simulations. Although the qualita- 
tive trend is clear, quantitatively it is still hard to com- 
pare the works of different lattice groups even for the 
same theory due to e.g. difference in the actions used. 
Only the recent data from Nf = to Nf = 12 from the 
same group [5] (shown as blue boxes in Fig[l]) allow for di- 
rect comparison. (These results are for chiral restoration, 
which for the fundamental quarks is believed to trace 
the deconfinement rather closely.) The main observation 
from Fig[T] is that the critical coupling j3 c changes by a 
substantial factor when the flavor number increases from 
Nf = to Nf = 12. Similarly, the N c = 2 theory with 
two adjoint fermions flows into the coupling at the in- 
frared fixed point l/g 2 = 0.20(4)(3) [13] . about a factor 
4 stronger than the critical one at the deconfinement of 
the N c = 2 pure gauge theory. 

2. Let us now turn to the mechanism of confinement, 
't Hooft and Mandelstam suggested a "dual superconduc- 
tor" model Q~7| related it to the Bose-Einstein Condensa- 
tion (BEC) of chromo-magnetic monopoles. For reviews 
of those ideas at T = in lattice context see e.g. [T]. 
Only recently it was realized that if monopoles are in- 
deed the emergent excitations one should better study 
them in a "normal" phase. Furthermore, in the region 
right above the transition T > T c such monopoles should 
be the dominant thermal degrees of freedom "ready" for 
BEC [IS] [H] ■ Lattice and model studies have shown that 
the effective coupling of the "magnetic plasma" does run 
as 1/g, inversely to the electric [3U1[3T]. Consequences of 
the magnetic scenario help to understand other lattice re- 
sults [23J as well as heavy ion experiments [2H US] . This 
study is based on the scenario that confinement occurs as 
Bose condensation of monopoles that are the dominant 
physical degrees of freedom in the plasma near T c . 

How can the light fermions affect the monopoles? It 
is known that the light fermions can become attached to 
them. The so called "fermionic zero modes" are specific 
bound states in which positive kinetic energy of localiza- 
tion exactly cancels the magnetic-moment-field interac- 
tion. (Their existence and number Nm are required by 
the topological index theorems, and thus insensitive to 
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any perturbative monopole deformations.) While such 
states may still be bosons, they carry flavor indices due 
to the fermions and are "distinguishable" , thus not con- 
tributing to the BEC of "unoccupied" monopoles. 

These zero modes are known explicitly for 't Hooft- 
Polyakov monopoles |26j which are present in gauge the- 
ories with adjoint scalars, such as in M = 2 supersym- 
metric gauge theories. Since flipping the charge and the 
spin leads to the same Dirac equation, antiquarks also 
have the same zero modes. Furthermore, for each of 
these zero mode states, it can be either populated or 
not, so the number of totally available states grows ex- 
ponentially ~ 2 2NfNM . Spectroscopy of those states in 
the supersymmetric setting was developed in 1990 's, see 
e.g. [35]. "Magnetic supermultiplets" have been explic- 
itly checked for two famous conformal theories, the Af=4 
SYM and the Af=2 SQCD with N f = 4: in bot cases 
one finds exactly the same set of spins/multiplicities as 
that in the original electric one. (Since those theories are 
electric-magnetic self-dual, their coupling cannot run at 
all!) There is one zero mode for the fundamental while 
two for the adjoint (Dirac) fermions: so these monopole- 
single-fermion states case into spin and | objects, re- 
spectively. While in the static case zero mode states are 
degenerate with the pure monopole, it is not so for the 
dynamical lattice monopoles with non-static paths. 

3. The BEC criteria for an interacting boson ensemble 
was proposed by Feynman 50 years ago for the study of 
liquid 4 iJe [29j and was recently generalized by the anal- 
ysis of Cristoforetti and Shuryak [30 . In the finite-T de- 
scription with periodic paths there appear "fc-clusters" of 
bosons interchanging their initial (at Matsubara time 0) 
and final (at Matsubara time j3 = 1/T) positions. Those 
clusters can be depicted as "Feynman polygons" with 
k points. Their probability depends on two competing 
factors, suffering from suppression due to the extra ac- 
tion exp{—kS ex ) while benefiting from large combinatory 
number of fc-polygons. The balance point marks the on- 
set of condensation by divergence of the sum over the 
/c-clusters. This method has been used in [22] where it 
was shown that lattice monopoles do have the BEC tran- 
sition exactly at T c . See also other studies of macroscopic 
clusters (or "percolation") such as[Tni[2n]- According to 
Feynman, quantitative BEC condition is a universal crit- 
ical value of the extra action per particle, which in three 
spatial dimensions is given by 

S ex < S c ~ 1.655 (2) 

Upon fulfilling {2}, long sequences of "hopping" bosons 
will occur, creating a macroscopic "supercurrent" . Its 
validity for interacting systems is demonstrated in |30j . 

Let us now apply the above criteria to the monopole 
condensation in QCD-like theories. The minimal ex- 
change action S ex for two nearest-neighbor bosons that 
are separated by a typical distance d = n^ 1 / 3 (with n the 
number density) during the Matsubara time from r = 
to r = (3 = 1/T could be estimated as 

S ex = to* + <P - m*f3 + S v (3) 




FIG. 1: Dependence of the critical lattice coupling /3 C at 
scale T c versus the number of fundamental quark flavors Nf 
in QCD-like theories. Blue boxes are from [9] with near- 
coincident boxes being lattice data for the same Nf with 
different N T which demonstrate lattice spacing consistency. 
Red diamonds are from various other literature. The thick 
blue line is the fitting curve, extended as dashed blue line be- 
yond Nf = 12. The black/purple/red dashed curves on the 
right are lines for vanishing beta function at 1,2,3-loop levels. 

with an explicitly written kinetic term, containing an 
effective mass to*, and the implicit potential term Sv- 
When close to condensation, the monopoles are very 
dense, with typical spatial separation d comparable or 
smaller than the inverse temperature /?, therefore justi- 
fying a further approximation of the above expression: 

S ex sa \m*Td 2 + S v = I (V) + S v . The term 

Sv = J^ T V(r(r))dT is basically the ratio of the inter- 
monopole interaction potential and T, also known as clas- 
sical plasma coupling Tm ^< V > /T ~ Sv As shown 
in Fig. 3b of [IB], at high T where the "magnetic scaling" 
d ~ l/(g 2 T) and g mag neUc ~ works, this ratio T M 
does not depend on the coupling or T and is a constant 
5, while close to condensation it decreases to a value 
of about 1 as T — >• T c for Nf = 0. The following onset 
condition will then be applied for the monopole conden- 
sation in QCD-like theories for the rest of our analysis: 

(f)(^Y i <S^2(S c -S v ) (4) 

The constant S c is of order one and its precise value is 
not needed as long as its iV/-dependence is negligible. 

4. Now, how would the transition get affected by 
adding light fermions? As already pointed out, the 
fermions can be attached to some of the monopoles via 
zero modes and effectively reduce the number of identical 
monopoles. Consider a monopole with one flavor of light 
quark added: for each of its allowed zero modes there is a 
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probability for it to be occupied by a fermion or not. Let 
us assume the ratio of the probabilities (occupied/not- 
occupied) to be / (a kind of zero- mode fugacity), we then 
see that the overall probability for a monopole (with Nm 
number of zero modes for each fermion flavor) to stay 
as a "pure" monopole is simply 1/(1 + J) 2Nm (with the 
factor 2 accounting for both quark and anti-quark contri- 
butions for Dirac fermions). So effectively the available 
density for BEC condensation will be n/(l + f) 2N f N **. 
Combined with the BEC condition in Eq.Q we obtain 
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(1 + /) 4Ar / Ar A-f/3 

{n/T 3 ) 2 / 3 



< S c 



(5) 



This implies that with increasing Nf, the density n has to 
increase and mass m* to decrease correspondingly so as 
to reach the same condensation condition. This pushes 
the transition to stronger coupling, therefore explaining 
the TVj-dependence of the critical coupling in FigfTj 

To make a semi-quantitative estimate, we use the 
following magnetic scaling relations that connect the 
monopole mass and density with gauge coupling: 
m*/T ~ l/g and n 1 / 3 /T ~ g 2 [HI OH HQ- Com- 
bined with the above condition, we obtain the criti- 
cal gauge coupling for monopole condensation to be: 
g(Nf) = g (1 + f) 4N f N M/ 15 where g is the correspond- 
ing critical coupling for pure gauge case. This can be 
further converted to the lattice couping (3 = 2N c /g 2 : 



P c {N f ) = A)(l + /) 



-8N f N M /15 



(6) 



As a concrete example let us focus on the case of 
fundamental fermions (with Nm = 1). The observ- 
able we examine is the critical (lattice) gauge coupling, 
(3 C ee (3{T = T c ) = 2N c /g 2 (T c ) as a function of Nf, as 
shown in figJT| In particular we concentrate on the data 
for Nf = to Nf — 12 from the same group [9], shown 
in Figjl] as blue boxes. We've used the above Eq.Q 
to make a fit for these data (blue boxes) and obtained 
the optimal value / f=s 0.154: the fitting curve is shown 
as the thick blue line in Figjl] Our model formula in 
Eq.(j6| with one parameter nicely describes all the data 
points in [9 from Nf — to Nf — 12! The suppres- 
sion factor / of monopole-quark as compared with pure 
monopole may be understood as follows: the monopole- 
quark has color-electric charge, and in the near-T c plasma 
it was known from previous studies [53] that the elec- 
tric particles are heavier than the magnetic particles by 
roughly AM ~ 2T C , thus leading to a suppression factor 
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-am/t _ 0.135 that is fairly close to / w 0.154. 



For completeness we've also displayed in Figjl the lat- 
tice results from various other groups [5HH1 [n! : these 
are shown as red diamonds. Admittedly, there are un- 
certainties due to different lattice actions and ambigui- 
ties associated with the possible differences between chi- 
ral/confinement transitions, which shall all be sorted out 
in the future lattice simulations. Nevertheless the quali- 
tative trend of decreasing lattice coupling with Nf is well 
in line with the data in [5] and with our model formula. 



The proposed mechanism suggests an approximate 
"NfNM scaling" , e.g. that the effect of adding N c funda- 
mental fermions is about as large as adding one adjoint. 
Essentially for a given gauge group, the critical couplings 
shall fall on one single curve when plotted against the 
combination NfN^j with various fermion representations 
and flavor numbers. Such scaling can be readily tested. 

5. To further "probe" the fermions' effect on con- 
finement transition, one may turn on a quark chemical 
potential fi q and see how the critical coupling changes 
accordingly. With the presence of a small fj, q , there 
will be difference if a zero mode is occupied by a quark 
or an anti-quark. This effect can be incorporated into 
the present model by replacing (1 + f) 2 in Eq.(|6| by 
(1 + /e 2 )(l + fe~ z ) with z = fJ, q /T. We then obtained: 
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*N f N M f 
15 (l + /) ; 



z 2 + 6(z 4 ) (7) 



Pc{N fl z) 
p c (N f ,z = 0) 

for small chemical potential /i q << T. Present lattice 
simulations, though not capable of handing finite /i di- 
rectly due to the sign problem, are actually able to ex- 
tract such dependence (or Taylor coefficients in fX q /T 
expansion) |32j and test the above prediction. Similar 
estimate can be made for the dependence on the isospin 
chemical potential as well as the axial chemical potential. 

The existence of these monopole-quark states also im- 
plies a contribution from them to the thermal fluctu- 
ations of the conserved charges (e.g. baryon number) 
they carry, thus contributing to the quark number sus- 
ceptibilities [331 [31]. A monopole may have one, two, 
or more zero modes occupied, but since the extracted 
/ « 0.154 is small the main contribution would be 
from the monopoles with one single zero-mode quark. 
The net baryonic density from these monopole-quark 
states (with the presence of a small quark chemical po- 



T 3 ■ 



2/ 



/T 3 



tential z = fi q /T) can be estimated as 
(n/T 3 ) (e z -e- z )f/[{l+fe z )(l+fe- z )] N f NM . This yields 
a contribution to the quark number susceptibilities Xi 
d^{n m _ q /T 3 )/d{ Pq /Ty^ as X r 9 



l+f) N f N > 



0.4—0.8 where we used lattice results for a total monopole 
density near T c to be n/T 3 - 2 - 4, and N f N M = 3. 
The so-obtained number makes a significant fraction of 
the lattice results around T c [33] ■ Higher-order suscepti- 
bilities can be estimated similarly and the multi-quark- 
monopole states with higher charges may be important 
there. Similar effects can be estimated along this line 
also for the isospin and electric charge fluctuations. 

6. It is well known that the N c = 2 theory is a 
very special case, with extra symmetry between quarks- 
antiquarks and mesons-diquarks. It also allows the finite 
density lattice simulations without the "sign problem". 
Lattice study of this theory was recently extended to the 
low-T finite-/^ region with Nf = 2,4 quarks by Hands, 
et al [35]. The quark density (per flavor) shown in their 
Fig. 3 displays a number of features: (i) a structure at 
|U R3 m^/2 as predicted by the rotation from the 'iptp to 
diquark condensate [3H1 [3Z] ; (ii) the usual quark Fermi 
sphere at higher //; and (iii) an unexpected growth of 
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quark density at still higher (i to about twice the value 
as expected from the Fermi sphere. The deconfinement 
as per the Polyakov loop appears concurrent with (iii). 

We now propose that this extra quark density in (iii) 
is due to the condensate of the monopole- quark states. 
The high quark chemical potential strongly favors states 
with quark numbers and efficiently converts the pure 
monopoles into monopole-quark states, thus explaining 
the deconfinement at about the same density. At such 
low T the dominant monopole-single-quark objects, be- 
ing bosonic, would appear mostly as a condensate, like 
the diquarks. Assuming standard effective potential with 
a repulsive binary interaction Vint = An 2 /4, one gets 
the condensate density growing linearly with fi, ubec ~ 
(// — m) /A at fj, > to, which is consistent with observa- 
tions of [35] . Furthermore, a similar density per flavor 
for both Nf — 2 and Nf = 4 is consistent with our 
view that these objects are dominantly states with only 
a single quark per monopole. States with multiple quark 
may become relevant at even higher density though. This 
proposal should and can be checked in many ways. Di- 
rect monopole-flavor correlations can be seen in the con- 
figurations of these simulations. One may also find the 



evaporation of this BEC as T is above certain critical 
value, into the Bose gas of such monopole-quark states. 

7. While we believe the main idea is robust, the 
model here is admittedly crude and intended to lead to- 
ward further studies and direct tests by dedicated lat- 
tice data. In particular, the correlation between the 
monopole line with flavor-carrying fermionic operators 
can be used to measure the probability for quarks "rid- 
ing" on monopoles. The discussion here is limited to the 
usual confinement, due to BEC of "empty" monopoles. 
More exotic objects, bosonic monopoles with quarks or 
even two monopoles bound by quarks, can in principle 
undergo BEC as well: see recent discussion in [55] . 

Note added: After submitting the paper we became 
aware of |39) discussing a similar problem in 2+1 dimen- 
sions. While the phase transition in this case is not of 
Bose condensation type and monopoles are substituted 
by vortices, the role of fermionic zero modes is similar 
and it also leads to reduced transition temperature. 
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